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At second order in perturbation theory, the r-modes of uniformly rotating stars include an ax- 
isymmetric part that can be identified with differential rotation of the background star. If one does 
not include radiation-reaction, the differential rotation is constant in time and has been computed 
by Sa. It has a gauge dependence associated with the family of time-independent perturbations that 
add differential rotation to the unperturbed equilibrium star: For stars with a barotropic equation 
of state, one can add to the time-independent second-order solution arbitrary differential rotation 
that is stratified on cylinders (that is a function of distance zu to the axis of rotation). We show 
here that the gravitational radiation-reaction force that drives the r-mode instability removes this 
gauge freedom: The exponentially growing differential rotation of the unstable second-order r-mode 
is unique. We derive a general expression for this rotation law for Newtonian models and evaluate 
it explicitly for slowly rotating models with polytropic equations of state. 


I. INTRODUCTION 

Unstable r-modes [li may limit the angular veloc¬ 
ity of old neutron stars spun up by accretion and may 
contribute to the spin-down of nascent neutron stars (see 
M for references and reviews). Spruit Q argued that 
angular momentum loss from the star would generate dif¬ 
ferential rotation, because the loss rate depends on the 
mode shape and varies over the star. Growing differ¬ 
ential rotation winds up and amplifies the star’s mag¬ 
netic field, and Rezzolla and collaborators [iSi , stud¬ 
ied the possibility that the energy lost to the magnetic 
field would damp out the r-mode instability. (In Spruit’s 
scenario, a buoyancy instability of the greatly enhanced 
magnetic field could power a y-ray burst.) To estimate 
the magnetic-field windup, Rezzolla et al. used a drift 
velocity of a fluid element; this is second-order in pertur¬ 
bation theory, but because the second-order velocity field 
had not been computed, they estimated it by integrating 
the first order velocity field. Subsequently, Cuofano et al. 
[Hill used this estimate of drift velocity to study the 
evolution of the r-mode instability damped by magnetic 
field wind-up.^ 

Following Spruit’s work. Levin and Ushomirsky found 
the differential rotation of the unstable r-mode in a toy 
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model of a spherical shell of fluid H- Sa [H then car¬ 
ried out the first computation of the differential rota¬ 
tion associated with a stable r-mode of uniformly ro¬ 
tating barotropic Newtonian stellar models and, with 
collaborators, looked at implications of the calculation 
for the unstable mode [3, [13 ■ The differential rotation 
arises at second order in perturbation theory as a time- 
independent, axisymmetric part of the solution to the 
perturbed Euler equations; for the r-mode whose linear 
part is associated with the angular harmonic Sa’s 
solution has the form 

= a^nCn (I)' (1) 

Here a measures the amplitude of the first-order pertur¬ 
bation, Cn is dimensionless and of order unity, the z-axis 
is the axis of rotation, and m is the distance from the axis. 
The function S),f'il{zu) is arbitrary. This ambiguity in the 
rotation law is present for the following reason. One can 
perturb a uniformly rotating barotropic star by adding 
differential rotation, changing the angular velocity from 
O to O -I- (5H(tJ7). If 5Ul{zj) is chosen to be quadratic in 
a, SUl{'Uj) = fl{'cu), it and the corresponding time- 

independent perturbations of density, pressure, and grav¬ 
itational potential $, constitute a solution to the time- 
independent second-order perturbation equations. Cao 
et al. [13 use a particular choice of to recompute 

the magnetic damping. 

In the present paper, we show that the second-order 
radiation-reaction force removes the ambiguity in the 
differential rotation associated with the Newtonian r- 
modes. In effect, the degeneracy in the space of zero- 
frequency solutions is broken by the radiation-reaction 
force, which picks out a unique differential rotation law 
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that depends on the neutron-star equation of state. We 
find an explicit formula for that rotation law for the un¬ 
stable r-modes of slowly rotating stars. 

To lowest nonvanishing post-Newtonian order, the 
growth time r of the radiation-reaction driven (CFS) in¬ 
stability of an r-mode is given by 


/3 = 


1 

r 


G 




where Cp is a dimensionless constant that depends on the 
equation of state. In using the Newtonian Euler equa¬ 
tion together with the radiation-reaction force at lowest 
nonvanishing post-Newtonian order, we are neglecting 
radiation-reaction terms smaller by factors of 0{R^l/c) 
and 0{GM/Rc^); this means, in particular, that we keep 
only terms linear in the dimensionless parameter (3/^. 

Three small parameters appear in the paper: The am¬ 
plitude a of the perturbation, the dimensionless growth 
rate and, in the final, slow-rotation part of the pa¬ 
per, the angular velocity il. For the logic of the paper, 
it is helpful to note that these three parameters can be 
regarded as independent of one another. The growth 
rate (3 can be varied by changing the equation of state of 
the material while keeping a and fixed; for example, 
in polytropes (stars based on the polytropic equation of 
state p = Kp^), one can change (3 by changing the poly¬ 
tropic constant K. 

The plan of the paper is as follows. Sect. m lists the 
equations governing a Newtonian star acted on by a post- 
Newtonian radiation-reaction force, with the star mod¬ 
eled as a self-gravitating perfect fluid. In Sect. nil we 
discuss first- and second-order perturbations of a uni¬ 
formly rotating star. From the second-order equations, 
we obtain a formal expression for the unique differen¬ 
tial rotation law of an unstable r-mode in terms of the 
first-order perturbations and second-order contributions 
that will turn out to be of higher-order in 17. Up to 
this point in the paper, the analysis holds for rapidly 
rotating stars. In Sect. IIVI we specialize to a slowly ro¬ 
tating background, keeping terms of lowest nonvanishing 
order in 17 and thereby obtaining an explicit formula for 
the radiation-reaction induced differential rotation. Fi¬ 
nally, a discussion section briefly comments on the valid¬ 
ity of the results for an accreting neutron star, when one 
includes magnetic fields, nonzero initial data for other 
modes, and viscosity. 

Our notation for fluid perturbations is chosen to make 
explicit the orders of the expansions in the amplitude a 
and angular velocity 17. The notation is defined as it is 
introduced in Secs. II and III, but, for easy reference, 
we also provide a table that summarizes the notation in 
Appendix|Al We use gravitational units, setting G = c = 
1 . 


II. NEWTONIAN STELLAR MODELS 


Let Q = {p, $} denote the collection of fields 

that determine the state of the fluid in a self-gravitating 
Newtonian stellar model. The quantity p represents the 
mass density, the fluid velocity, p the pressure, and <I> 
the gravitational potential. For a barotropic equation of 
state p = p{p), the specific enthalpy h of the fluid is 


h = 


Jo P ’ 


and we define a potential U by 


U = h + ^. 


( 2 ) 

( 3 ) 


The evolution of the fluid is determined by Euler’s equa¬ 
tion, the mass-conservation law, and the Poisson equa¬ 
tion for the Newtonian gravitational potential. These 
equations may be written as 

+ + = (4) 

0 = dtp + Va{pv^), (5) 

= 47rp. (6) 

The version of the Euler equation that we use, 
Eq. ([1]), includes fcR, the post-Newtonian gravitational 
radiation-reaction force (per unit mass). This force plays 
a central role in the nonlinear evolution of the r-modes 
that is the primary focus of our paper. It is given by 

- ^ ^ (-1)^+1^ 

Jgr-2^ 2^ ^ -/I 

Z>2 |m|<i I 

2r^Y^rn ^21+2gem 2v X V(r^F^'") \ 

~^TTT Vl 


where 5R(F) denotes the real part of a complex quan¬ 
tity Z. The quantities and are the complex 
mass and current multiple moments of the fluid source 
(cf. Thorne [l^ Eqs. 5.18a,b) defined by, 


j£.m 

^im 


VI 


pr Y 




2m 

VITT 


pVv ■ 




with Ni the constant 


16^ l{i + 2){£+l) 
{2i + iy.\\l 2{e-i) 


( 8 ) 

(9) 


( 10 ) 


The functions Y^"^ are the standard spherical harmonics, 
while the FJ™ are the magnetic-type vector harmonics 
defined by 


-^£m 
^ B 


r X VF'™ 

VW+^' 


( 11 ) 
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We use the normalizations 1 = J \ Y^"^\‘^dcos 9dcj) and 
1 = J \ Y^"^\^d cos 9d(f> for these spherical harmonics. 
In Cartesian coordinates f is given by r = {x,y,z). 
We point out that this expression for the gravitational 
radiation-reaction force, Eq. agrees with the mass- 
multipole part of the force given by Ipser and Lind- 
blom [ 2 ^. It also agrees with the current-multipole part 
of the force given by Lindblom, et al. (following 
Blanchet [ 2 ^ and Rezzolla, et al. [l^) for the £ = 2 
and m = 2 case. The general form of the force given in 
Eq. ([71), however, is new. 

The post-Newtonian radiation-reaction force is gauge 
dependent, so the expression for it is not unique. We 
derived the expression for the force given in Eq. © by 
requiring that it implies a time-averaged (over several 
oscillation periods) power {{dE/dt))\GR (which is gauge 
invariant), and angular momentum flux {{dJ/dt))\GR lost 
to gravitational waves that agree with the standard post- 
Newtonian expressions, cf. Thorne [l^. We present ex¬ 
pressions for these flux quantities in Appendix|B]that are 
equivalent to, but are somewhat simpler than the stan¬ 
dard ones. 

We consider small perturbations of rigidly rotating, ax- 
isymmetric, barotropic equilibrium models (models with 
a barotropic equation of state). The fluid velocity in 
these equilibria is denoted 

V = £1$^ (12) 


These derivatives can be used to define a formal power 
series expansion for 6Q\ 

SQ{a,t,x) = aS^^'^Q(t,x) + a‘^ d^‘^'^Q{t,x) + 0{a^). (16) 


Each point in the interior of the unperturbed star is, for 
sufficiently small a, in the interior of the perturbed star; 
the derivatives defined in Eq. (USD and the formal 
power series expansion in Eq. dUD are thus well-defined 
at all points of the interior of the unperturbed star, but 
may diverge at the surface. We consider constant-mass 
sequences of stellar models, i.e., models whose exact mass 
perturbations, SM = M{a) — M{a = 0) vanish iden¬ 
tically for all values of a. The integrals of the n**'- 
order density perturbations therefore vanish identically 
for these models: 


1 d”M(a) 

nl da" 


Q;=0 


/ 


S^^^p^d^x. 


(17) 


The exact (to all orders in the perturbation parameter 
a) perturbed evolution equations for these stellar models 
can be written in the form 


SE‘^ = (5t-h fl£0)du“-H -H V“(5C/, 

+ 6v^Wi,6v‘^ = (18) 

0 = {dt^£.(p)Sp^a{pSp Sv^"), ( 19 ) 

= AnSp, ( 20 ) 


where (p generates rotations about the z axis; in Cartesian 
coordinates, 4> = (—y,x, 0). Eor barotropic equilibria, 
Euler’s equation reduces to 

0 = Wa{h + ^- (13) 

where h is the specific enthalpy of the fluid and w is the 
cylindrical radial coordinate, = x"^ + . The surface 

of the star is the boundary where the pressure and the 
enthalpy vanish: p = h = 0. 


III. PERTURBED STELLAR MODELS 


We denote by Q{a, t, x) a one-parameter family of stel¬ 
lar models. Eor each value of the parameter a, Q{a, t, x) 
satisfies the full nonlinear time-dependent Eqs. ©-©■ 
We assume that the model with a = 0 is an axisymmetric 
equilibrium model, as described in Eqs. dUll-dlSD. The 
exact perturbation 5Q, defined as the difference between 
Q{a) and <5(0), is defined everywhere on the intersection 
of the domains where Q{a) and <5(0) are defined: 


5Q{a, t, x) = Q(a, t, x) — Q{0, t, x). (14) 

It is also be useful to define the derivatives of the 

one parameter family Q{a) evaluated at the unperturbed 
stellar model, where a = 0: 


5^^^Q{t,x) 


1 Q{a, t, x) 
n\ 9a" 


(15) 


where is the Lie derivative along the vector field 0, 
and p is the density of the unperturbed star. The ex¬ 
act perturbed gravitational radiation-reaction force 5Pgr 
that appears in Eq. dUD is given by 


0JGR = 2_^ 2_^ -r:::- 


327r 






l>2 \m\<l 


1 dt'^^+'^ 

2Sv X V(/y^'") I 


VI 




where 


gjlm = 5p/Y*^'^d^X, 




2Nt 

VITT 


pSv + dp + sVj 


( 21 ) 

( 22 ) 




(23) 


It is convenient to decompose the perturbations 5Q 
into parts 5mQ that satisfy the pure Newtonian evolu¬ 
tion equations, and parts 5b.Q caused by the addition of 
the gravitational radiation-reaction force. In particular 
the non-radiative stellar perturbations 5nQ satisfy the 
perturbed Euler equation: 

5E = 0. 


01=0 


(24) 
















4 


When the effects of gravitational radiation-reaction are 
included, the complete perturbation, SQ, satisfies the 
Euler equation driven by the gravitational radiation- 
reaction force 


SE = SfcR- 


(25) 


A. First Order Perturbations 


The classical first-order (in powers of a) r-modes have 
angular and temporal dependence sm 

= ^N^P- sintjjN, (26) 

sinV'AT-h cosV’AT, 

(27) 

4'^^^ = 4^^^- sinV-iv, (28) 

= 4^^- simjjN, (29) 

where ipN = toNt + m(j), with m ^ 0. The tensor 

P\ = 5\-vj-^r(l>b (30) 

is the projection operator orthogonal to ()>“, and = 

S^^'^Q{va, z) depends on the cylindrical coordinates zu and 
z, but not on ^ or t. The origin of time has been chosen 
to give the perturbations definite parity under the diffeo- 
morphism (j) —>■ —cj) at t = 0. We use the term (j)-parity to 
mean parity under this transformation. The subscripts ± 
indicate that and 4^^- parts of odd 

0-parity scalars, while is part of an even 0-parity 

vector field. 

When gravitational radiation reaction is included, the 
Euler equation is altered by the relatively weak radiation- 
reaction force fcR- The first order radiation-reaction 
force can be written in the form: 


changed in the following way (schematically): 




{ 5 ^^^P- sin'0 -I- 4^/5+ cos'0 e^*, 

(32) 


4^^- 


w 0“0b COS0 -I- P“b sin0 




(4^4+ 4^4) X 




j(i)$ 


G7 0“0{, sin0-I-P“f, cosV' 
^4^^- +4^^-) sin0e^* 




^4^^- + 4^^-) sin0e' 


pt 


-l-4^®-i- costpe^*, 


pt 


(33) 

(34) 

(35) 


where tp = 'ipN + tpn = ujnI + ujRt + m0. The radia¬ 
tive corrections are smaller than the non-radiative 

perturbations 44 by terms of order 0{/3/ujn). The ra¬ 
diative correction ujn to the frequency, is smaller than lon 
by a term of order 0{/3/ujn)^, so we ignore that change 
here, setting ip = ipN-'^ 


The radiative corrections to the r-mode, 44i are de¬ 
termined by substituting Eqs. dsa)-® into the first- 
order perturbed mass conservation and Euler equa¬ 
tions. After applying the equations satisfied by the non- 
radiative parts of the perturbations, 44 j the resulting 
system of equations can be divided into parts propor¬ 
tional to sin0jv and cos "0 at respectively, each of which 
must vanish separately. The resulting equations can be 
divided further into a set that determines 4^^/5-j 44-, 
and 4^4’ another that determines 5^^p+, 44+, 
and 4^^-- 


= 044 + + 4VGiJ+, (31) 

where 0 is the growth rate of the r-mode instability, 
and 4^/Gfl-i- is (by definition) the even 0-parity part of 

the radiation-reaction force that is orthogonal to 4^+ 
and that therefore does not contribute directly to the 
energy evolution of the mode. Equation m implies 
that the odd 0-parity part of the radiation-reaction force, 
fcR-j vanishes when the classical r-mode is chosen 
to have the 0-parity given in Eqs. (|26|l -([29 |) . The gravi¬ 
tational radiation-reaction force causes an instability by 
introducing an imaginary part 0 to the frequency of the 
mode. The overall structure of the modes is therefore 


The equations that determine the radiative corrections 
having the same 0-parity as the classical non-radiative r- 


^ Friedman and Schutz fSsj derive the following general expression 
for the frequencies of the modes of Lagrangian systems (including 
Newtonian fluids with gravitational radiation-reaction forces): 
0 = A{u -(- i/3)^ — (5 -h iD){u) -h ip) — C, where A, S, C and 
D are real. The term D vanishes for non-dissipative Newtonian 
fluid stars. When D is small, it is straightforword to show that 
the real part of the frequency, oj, differs from the frequency of 
the non-dissipative D = 0 system, by terms of order D^\ 
u) = u)^ -h O(D^). It is also easy to show that the imaginary 
part of the frequency 0 is proportional to D for a mode with 
= 0. 
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modes are then 


(wjv + mO) + mptzj (pa v[ 


(!)„' 


+V, 


, (pP%(5Wz)t) = 0, 


(36) 


[(a;AT + mQ.)(j)a + 2wQ.Va'^] = —rnd}^ U-^ (37) 

2 

{ujn + mrt)P^b H- Q.V°^w(j)h 

VD 




= P“'’Vb4^’[7_. 


(38) 


These equations are homogeneous and are identical to 
those satisfied by the classical r-modes. The solutions 
for and are therefore proportional 

to the classical r-modes: 5^^/3_, and . The 

effect of adding these radiative corrections to the classical 
r-modes is simply to re-scale its amplitude. We choose to 
keep the amplitude, a, of the mode fixed, and therefore 
without loss of generality we set 

0 = 5^n^P- = S^^'>U-=S^^'>vl. (39) 


It follows that the first-order radiative corrections have 
((eparity opposite to that of the classical r-modes: = 

(5^^P+, = 5^^17+, and ■ They are 

determined by the equations 


{uj]\f + mil) p + mptu 


(1),- 


-V, 


(pP“f,4){)') =/35^V, (40) 


[(wat -I- mVt)pa — 2wilVa'^] + W ^''rU 








{ujn + m^)P^b - 

w 


(41) 


(42) 


The general solution to the inhomogeneous system, 
Eqs. (Hn]) - (H^ . for (5j^^p, and consists of 

an arbitrary solution to the homogeneous equations (ob¬ 
tained by setting p = <5+Vga = 0 ) plus a particular 
solution. These homogeneous equations are identical to 
Eqs. (1551) - (1551) . so their general solution is a multiple of 
the classical r-modes. Because their ^-parity is opposite 
to that of the classical r-modes the effect of the homoge¬ 
neous contributions ^^rU, and 4^’ is to change 
the overall phase of the mode. We choose (by appropri¬ 
ately adjusting the time that we label < = 0) to keep this 
phase unchanged, and we can therefore, without loss of 
generality, set to zero the homogeneous parts of the so¬ 
lutions to Eqs. dSSD-dMl). The inhomogeneous terms on 
the right sides of Eqs. (15ni) - (H51) . p and 5^^ Jqr^ are 
all of order /3. Thus the particular solution to Eqs. (gO])- 
(155)) must also be of order /3 as well. It follows that the 
radiation-reaction corrections to the first-order r-modes 
4^(3 are smaller than the classical r-modes 4^Q by 


terms of order 0[l3/w). To lowest-order in /3, therefore, 
the corrections to the first-order r-modes in Eqs. (1551) - 
(l35)) simply change the overall scale of the mode by the 
factor e^‘: 


B. Second-Order Perturbations 


The second-order perturbation equations are a sum 
of terms linear in and terms quadratic in 

For example, the second-order perturbation of the Eu- 

1 (f 

ler equation, = --r-^E° 


2 da '^" 


includes the term 


a=0 


which serves as an effective source term 
for the second-order perturbations and In 

the absence of gravitational radiation reaction, it follows 


( 2 ) 

that the second-order Newtonian r-mode 6 )^ Q is a sum 
of terms of three kinds: a term with angular and tem¬ 
poral dependence cos(2-0Ar), where tpN = 'mfp + ^Art, a 
term with dependence sin(2'0Af), and a term that is time 
independent and axisymmetric. This time-independent 
axisymmetric part of the velocity perturbation can be 
regarded as differential rotation. As we have emphasized 
in the Introduction, the second-order Newtonian r-modes 
are not determined uniquely: Given a particular solu- 
tion 5)^pQ to the second-order Newtonian perturbation 
equations with perturbed velocity field 44'*; there is 


f 2) 

a family of solutions 5p; Q with perturbed velocity field 


4 ^^^“ = 5^^'^pV°' + , where 4^^(^) i® arbi¬ 

trary. This degeneracy is broken by gravitational radi¬ 
ation reaction. The presence of the radiation-reaction 
force picks out a unique that displays the gravi¬ 

tational radiation driven growth of the second-order r- 
modes: oc 

To find this differential rotation law, one must solve the 
second-order axisymmetric perturbation equations with 
radiation-reaction force for the axisymmetric parts of the 
second-order r-modes. Denote the axisymmetric part of 
a perturbation 6 Q by (SQ), and denote by the ex¬ 

ponentially growing differential rotation of the unstable 
r-mode: 


= [(444 + 4^f^(^)]e"^*- (43) 


Without solving the full system, however, one can ob¬ 
tain a formal expression for in terms of the known 

first-order perturbation together with other parts of the 
second-order axisymmetric perturbation. As we will see 
in the next section, this expression is all that is needed to 
find to lowest nonvanishing order in O: The other 

parts of the second-order perturbation give only higher- 
order contributions. Finding this formal expression for 
and showing that it is unique are the goals of the 
present section. 

We now turn our attention to solving the perturbation 
equations for the axisymmetric parts of the second-order 
r-modes. The axisymmetric parts of the second-order 
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perturbations can be written in terms of their radiative 
and non-radiative pieces: 


The non-radiative parts of the perturbations 

are determined, up to a perturbation that adds differen- 


(^(^V) = 

((4V) + (4V)) 

e^3\ 

(44a) 

( 2 ") 

tial rotation by the axisymmetric parts of the 

Newtonian Euler and mass-conservation equations: 

( 5 ( 2 ) 1 ;“) = 

((4^^“) + (4^^“ 


(44b) 

+ v“(544 = -(54 


(5(2)C/) = 

((544 + ( 5(^)4 


(44c) 



(52) 

( 5 ( 2 )$) = 

((4^4+ (4^4 


(44d) 

Va 

p(54u“) + (54p54^;“)' 

= 0. (53) 

{s^^^roR) = 

{s^R^foRy^'- 


(44e) 

( 2 ^ 

Given a particular solution 6 }^pQ to these equations, we 


These quantities are determined by the second-order ax- 
isymmetric parts of the perturbed stellar evolution equa¬ 
tions: 

2/3(,5(2)u“) -b -b 

= (45) 

2p{S^^'>p) + U(5(2)u“) -b = 0, (46) 

=47r(5(2)p). (47) 

The uniqueness of the second-order differential rota¬ 
tion can be seen as follows. Let and 

be two solutions to the second-order perturba¬ 
tion equations, Eqs. (|45|), dUl), and gzl), associated with 
the same time-dependence and with the same first- 

order solution The difference {S^^'>Q) — {S^^'>Q) 

of the two solutions then satisfies the linearized Poisson 
equation and the linearized Euler and mass conservation 
equations obtained by setting to zero the terms involv¬ 
ing and Eqs. (HSl) and (HSl) . That is, 

((5(2)q) _ (5(2)Q))e2/3i 

is an axisymmetric solution to 
the first-order Newtonian perturbation equations. But 
the Newtonian star has no such solution, no mode with 
growth rate 2/3. Thus {{S^^'>Q) — {S^'^^Q))e'^^* = 0, im¬ 
plying that is unique. (Note, however, that the 

decomposition (1431) is not unique: The arbitrariness in 
the differential rotation of the Newtonian r-mode means 
that one is free to add to an arbitrary func- 

tion /(to) if one simultaneously changes to 

( 2 ") ( 2 ) 

We now obtain equations for Q and S)^ Q. Keeping 
terms to first order in /3, the terms quadratic in first-order 
perturbed quantities that appear in Eqs. dlSl) and (gS) 
have the forms. 




On ^ V bON ^ / 


+13 {Sr 


( 2 ) 






-'N 




where 

/3(4^V“> = + 


(48) 

p2/3t 
-• 1 

(49) 

>, (50) 
(51) 


-'N 


differential rotation of Eq. (1431) that is picked out by the 
radiation-reaction. 

We define the radiative part of the perturbation, 
(44) j by requiring that it be created entirely by the ra¬ 
diation reaction forces; (44) is therefore proportional 

f2l 

to the radiation reaction rate /3. When {S}^'Q) satisfies 
the Newtonian equations (l52l) and (l53)) , the axisymmetric 
parts of the full perturbed Euler and mass-conservation 
equations with radiation-reaction have at 0(/3) the form 

2/3(4^^'“) + 211(4^-b V“(4V) 

= ( 44 “ 4 -/?( 44 “), ( 54 ) 

Va(p(4^«“)) =-2/3(4V)-/3Va(4V“). (55) 

To find an expression for S}^ we first write 

(4^'*^“) ^'S (44”^) + 44(to)(/)“ and move the term 
involving (44°) i-be right side of Eq. (1551) : 


(56) 


2/344(to)(/)“- b 2fl(4^w'’)Vb(/)“ -b V“(44) 

= /3(4^E“), 

where 

^(4)f“) = (4Vgh) - 2/3(44°) - 13 (44“). (57) 

We next write the components of the axisymmet¬ 
ric part of the second-order perturbed Euler equation, 
Eq. (1561) . in cylindrical coordinates: 

2/3to44(to) -b 211(4’?^'^) =/3 w( 44'^), (58a) 
-2flTO(4^^^'^) = -5.^(44)-b/3(4^E'^), (58b) 

0 = _a,(4)c/)+/3(4)F^). (58c) 

Using Eq. (I58al) to determine (4^^^)i tbe axisymmetric 
part of the second-order mass conservation Eq. dSS]) can 
be written as 


/3 


2Uto 

+5, 




pw^ (^(4^^“^) — 24^^(^) 

p{SrV^)\ = -2^(4V) -PVa{5fw^). 

(59) 
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The star’s surface is defined as the p = 0 surface. Be¬ 
cause 5(2) 

p is a derivative evaluated at a = 0, it has 
support on the unperturbed star. While the density per¬ 
turbation p is not finite for some equations of state at 
the surface of the star, it is integrable in the sense that 
6 ^^'^ f pdz is finite, as one would expect from the inte- 
grability of the mass-conservation condition in Eq. (1171) . 
In particular, for polytropes with fractional polytropic 
index 0 < n < 2, p diverges at z = zs, but, as we 
show in Appendix O 6 ^^^ f pdz is finite. Here we de¬ 
note by zsivu) the value of z (the Cartesian coordinate 
axis parallel to the rotation axis) at the surface of the 
unperturbed star. 

We now multiply the second-order mass conservation 
equation, Eq. dMI), by 2 wflj(3 and integrate with respect 
to z over the support of the star. It will be convenient to 
extend the domain of integration to extend slightly be¬ 
yond the surface of the unperturbed star. Because each 
integrand has support on the unperturbed star, we sim¬ 
ply take the integrals to extend from — oo to oo instead 
of —Zs to zs- We then have 


0 = 


/ OO 

dz(5^^'^ p'^ 

-OO 

/ OO 

-OO ^ 

/ oo 

-OO 


(60) 


The second integral on the right side of Eq. (l60l) can be 
re-written as 


dz9„ 


pvj^ 


/ OO 

dzpw'^ - 2(5^^H(ro)j . (61) 


of area ■wd(l)dz. The term is then given by 



/ OO 

-OO 





dz(4^V“). 


(63) 


With this simplification, Eq. (1621) can be written in the 
form: 


2n7^d^^H(cc7) 



dz p = 



dzp{5fF^) 





dw'z 



dz{5%'^p) 


-|-2tJ7fl 



dz(4^%^). 


(64) 


This provides a formal expression for in 

terms of the first-order perturbations that comprise 
and and the second-order perturba¬ 

tion (s^np)- ^ 

Together with it determines the differential 

rotation of the unstable r-mode. 

We conclude this section with a discussion of two sim- 
plifications in evaluating ^{w), one from the fact that 
we work to first order in the growth rate /3, the second 
from the slow-rotation approximation of the next section. 
The first is a simplification of the expression Eq. (l64l) for 
the radiation-reaction force. The integrand of the first 
term in Eq. is given by the (()-component 

of Eq. dlZl): 

= (5^%^^) - 2f3{S^^\*) - (3 {Sfv^). (65) 


The expression in Eq. (1601) can then be integrated from 
m = 0 to vu, using Eq. m, to obtain an expression for 

/ OO pOO 

dz p = zu'^ / dzpl^S^^F^^ 

-OO J —oo 

pvs poo 

-|-4f2 / dw'w' / dz(5^^'^p') 

Jo J — oo 

pVJ poo 

+2n / dvu'w' / dzVa( 4 ^’lE“). (62) 

Jo J — oo 

Because of the axisymmetry of its integrand, the third 
term on the right side of Eq. (15^ is, up to a factor of 
27 r, the volume integral of a divergence. The boundary 
of the three-dimensional region of integration has two 
parts: One is just outside the surface of the star, where 
vanishes; the second is the cylinder at constant vo 
from —Zs to Z 5 , with outward normal VoCiJ and element 


To evaluate (<^pVgp)i must find the axisymmetric, 
second-order, part of the expression for 6 fan on the 
right side of Eq. (ED. Recall that the axisymmetric 
parts of any second-order quantity have time dependence 
g2/3t^ The first three terms in the bracketed expression in 
Eq. (|2T]) involve high-order time derivatives of or 

i5(2)S'^°, and are therefore proportional to high powers of 
P and can be neglected. We are left with only the fourth 


® As mentioned above, Appendix [Clshows that assuming smooth¬ 
ness of the displacement of the surface as a function of a and x 
( 2 '\ 

implies integrability of A simpler way to see that the right 

side of Eq. (63) is finite is to note that smoothness of the displace¬ 
ment of the surface implies one-sided differentiability of at 

the surface. The perturbed mass conservation equation, Eq. dsg, 
then implies that the combination 2{5^^p} H- is fi¬ 

nite at the surface and hence integrable. This is enough to imply 
that the expression in Eq. (I64II for is finite. 







term, 


{S^rIgr) 


{-lYm 

8 tt\/£ 


x5r(4^Vx v(/r") 


d2e+igWsU 




( 66 ) 


The second simplification involves the quantities 
and {S^rW^) that appear in Eq. (IMl) . They are 
defined in Eqs. (1501) and m- Using the general expres¬ 
sions for the first order perturbations given in Eqs. (I32[) - 
(I35D , we can express these quantities in terms of the first 
order perturbations: 


(/34^V“) 





(67) 


( 68 ) 


As we will see in the following section, these terms and 
the term involving p in Eq. (1641) are higher order in 
than the first two terms of Eq. (1551) and can therefore be 
neglected when evaluating for slowly rotating 

stars using Eq. (l64l) . This fact is essential, because 
itself depends on 4^^- 

This discussion has been somewhat abstract but quite 
general. Apart from assuming the integrability of the 
perturbed density so that mass conservation, Eq. (II3, 
can be enforced, no assumption has been made up to 
this point about the particular equation of state of the 
matter in these stellar models, nor has any assumption 
been made about the magnitude of the angular velocity 
of the star. In order to proceed further, however, we will 
need to assume that the stellar model is slowly rotat¬ 
ing in a suitable sense. To find an explicit solution for 
dY'^D,(w), we will also need to make some choice for the 
equation of state for the stellar matter. The slow rotation 
assumption and its implications are discussed in Sec. El 
while the complete solution for the second-order 

r-mode angular velocity that is driven by gravitational 
radiation reaction, is determined in Sec. lYl for the case 
stars of composed of matter with a range of polytropic 
equations of state. 


IV. SLOW ROTATION EXPANSION 


We consider the one-parameter families of stars Q = 
(5(0) composed of matter with a fixed equation of state, 
and having masses that are independent of the angular 
velocity: M(0) = Mq. The structures of slowly rotating 
stellar models in these families are conveniently written 
as expansions in the dimensionless angular velocity, 


O = 


O 

On 


8 

where Oq = a/ and Mq is the mass and R the 
radius of the non-rotating star in the sequence. The slow 
rotation expansion of these stellar models is denoted, 

Q = ^g„0” = Qo + Qi^ + Q202 + 0(03). (70) 

n—0 

For equilibrium rotating stars these expansions of the 


basic fluid variables have the forms: 

p = po+P2 02 + 0(04), (71) 

= 0(^“, (72) 

P = Po d- P2 O^ -l- 0(0^), (73) 

$ = $o + «'20^ + 0(0'‘). (74) 

We will represent the perturbations of these stellar mod¬ 
els SQ as dual expansions in the mode amplitude a and 
the angular velocity parameter O: 

= (75) 

n,k 


Our main goal here is to determine to lowest-order in 
angular velocity the axisymmetric part of the second- 
order perturbations of the r-mode angular velocity field 
that is driven by the gravitational-radiation in¬ 
stability. Doing this requires the explicit slow-rotation 
forms of the first and the second-order perturbations. 
These slow-rotation expansions are described in the re¬ 
mainder of this section. 


A. First Order Perturbations 


The effect of the first-order gravitational radiation- 
reaction force on the structure of the classical 

r-mode (beyond its overall effect on its amplitude) was 
first studied (for i = 2) hy Dias and Sa [T^. We agree 
with the results they obtain but will need to clarify their 
meaning. We also extend the calculation to general val¬ 
ues of i. 

To first order in mode amplitude a and lowest non¬ 
trivial order in angular velocity D, the classical r-modes 
with the ^-parity described in Sec. IIII Al can be written 
the form 


j(i) 


N 7*1 

x(i) - 

oyvi 


= = 0 , 

^^0 (sin^ 


(76) 


i Vi?y 


(77) 


where ^{Z) is the imaginary part of a quantity Z. An 
equivalent expression for the classical r-mode velocity in 
terms of vector spherical harmonics is 


On vi 


3 , 


iAer^ 

2 


\/-££ Aut 

1B e 


(-l)^Y^-^e-*‘^‘ 


(78) 

,(79) 


(69) 
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where Ag, is given by 


A, = (-l/2'(£-l)!. 


UTTi{£+l) .^1 

( 2 ^+ 1 )! 


^no. (80) 


The frequencies of these classical r-modes have the form 

{£-!)(£+ 2)^ , 


UJn = 


1 


+ (81) 


At this order in il, the r-modes do not affect the fluid 
variables Sp and 6p, which are Because of this, 

the r-mode velocity field at order does not depend on 
the equation of state. 

Four features of the gravitational radiation-reaction 
force are important in determining the way it alters each 
r-mode: a) The ())-parity of fcR, as shown in the last 
section, is opposite to that of the classical mode; b) its 
magnitude, as shown below, is dominated by the current 
current multipole c) it can be decomposed in the 
manner 




(82) 


where the two terms in the decomposition are orthog¬ 
onal with respect to a density-weighted inner product, 
/ y/g X pq 6^^'^V ■ fcR = 0; and d) as we show be¬ 
low, 5^1^ fcR is a gradient, fcR = 

It is straightforward to evaluate the multipole moments 
of the r-modes using Eqs. (l2^ and (l23l) and the expres¬ 
sions for the classical r-modes from Eqs. (dl) and (1771). 
The expressions for the non-vanishing multipole moments 
of the r-modes can be written in the form 




= (-1)^(5).)^ 5 


.AgNge’-' 


N 

iujt pR 


VITT 


,.2i+2 


Po dr. 


(83) 


Inserting these expressions into the formula for the grav¬ 
itational radiation-reaction force, Eq. (ED), we find 


foR = 


{-l)^Ng 

8 tt 


5R 


[vr+T ^ 


n 


+ £A-4> X V(rV") 

VI 




(84) 


This expression can be re-written as a linear combination 
of VyV V(r^YV using the identity 

$ X V{VyV = iVW+TiVYV - zV{r^YV- (85) 

The resulting expression for foR can therefore be 
written in the following way: 


S^^^Igr = 


where j3 is given by 


rR 


J21+2 


A^{£^-!){£ +2) 


podr, 


( 86 ) 


(87) 


and where fa r is defined by 


S^l^fGR = - 




.2^+2 


Stt 

£ + 1 
?(i), 


pdr 


5? 




vW+T) 


( 88 ) 


This expression for fan can be rewritten as a gradi¬ 
ent. 


S^l^fGR = ^{iPAi 


=: 


[r^+i cos0r"e“‘] | 


(89) 


Eqs. (1861) and (1891) give the decomposition of Eq. 
and the orthogonality of the two parts. 


Pq 5^^V ■ 6^1"^ faR^gd^x = Q, (90) 


is implied by the relation 


^abc 


Va(cos 9 Y“)^brS/^Y^^^d^ x 


= - e 


/' 


ibc , 


cos6»F"Vf,rVaVcF“v^d^x = 0, 


(91) 

where y/gd'^x is the volume element on the sphere: 
y/gd'^x = — r'^dcos 9 d(j). At this order in fl, the den¬ 
sity Po plays no role in the orthogonality, but it is with 
respect to the density-weighted inner product that the 
operators appearing in the perturbed Euler equation are 
formally self-adjoint. 

It follows that faR is the part of the gravitational 
radiation-reaction force that does not contribute directly 
to the exponential growth of the classical r-mode insta¬ 
bility and that the coefficient £3 is the growth rate of the 
gravitational radiation driven instability in the r-modes. 
Substituting into Eq. (1871) the expressions for Ng from 
Eq. (|I0p and the r-mode frequency uj from Eq. ED gives 

|(2f+l)!!p {e+ij f 

which agrees with the expression for the gravitational 
radiation growth rate of the r-mode instability given in 
Lindblom, Owen and Morsink Q. 

These expressions for the slow rotation limits of the 
radiation-reaction force confirm the general expressions, 
e.g. Eq. (1311) . used in our discussion of the general proper¬ 
ties of the first-order r-modes in Sec. lIII Al It follows from 
that discussion that the general form of the first-order r- 
mode velocity, to lowest order in the angular velocity of 
the star, is given by 


sVvi e^*. 


(93) 
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To evaluate fl using Eq. (l64l) , we need to determine 
and or at least to show that they are negligi¬ 

bly small compared to other terms in the equation. We 
show in the heuristic argument below that = 0(/3fi) 
and = 0(/3fi^), which will allow us to neglect them 
in our slow rotation expansion. A more precise version of 
the argument is given in Appendix ID] The fact that 
is higher-order in H than <5^^ is the reverse of their re¬ 
lation in the classical r-modes. This reversal depends on 
the appearance of the gradient in the decomposi¬ 

tion of the gravitational radiation-reaction force 

The equations that determine Eqs. (gni)-®, 

can be written more compactly as 

(ccjv + ^fl)4^p + V- (94) 

(ujjY + £ 

= (95) 

The value of is fixed by the curl of the perturbed 
Euler equation, (ESI): 


V X 


(wat + £+ 2Vt5''^’v ■ Vcj) 


(1): 


= 0 , 


(96) 


which involves only Its two independent compo¬ 

nents give two relations for the three components of 
in which all coefficients are 0(17). All components of 
are therefore of the same order in fl. Similarly, the 
two relations among and S^^p given by the 

equation of state and the Poisson equation imply that 
and are of the same order in fl. The continu¬ 
ity equation, (IMI) . then implies that = 0(fl(5^^p). 
Finally, the (/(-component of the Euler equation gives, to 
lowest order in fl, 

= <5Wjr -H C>(fl2<5Wp). (97) 


From its definition in Eq. (15^ it follows that = 

0{n/3), which then implies that = O^pVt) and 

5^^^v = om^) 

Dias and Sa [TJ find, for an £ = 2 perturbation, a so¬ 
lution that is a sum of a) our solution with 

5^rU given by Eq. (DTI) and b) a solution to the homo¬ 
geneous equations with (/(-parity opposite to that of the 
Newtonian r-mode As noted above, adding part 

b) of their solution is equivalent to changing the initial 
phase of the perturbation. 


B. Second Order Axisymmetric Perturbations 

In computing the quadratic terms that enter the 
second-order perturbation equations, it will be useful to 


have explicit expressions for the classical r-mode (J^^i 
in cylindrical coordinates (tn, z, (p), 

On 

( zu \ 

= -noz\^—j cos(£(/( + a;£). 

(98a) 

On Oi 

/ ZU \ ^ 

= floR\^—j cos{i(j) + cot), 

(98b) 

On Vi 

= sin(£(/(-b wt). 

(98c) 


From these one finds explicit expressions for the cylindri¬ 
cal components of the quadratic term 
which appears as a source in the second-order Euler equa¬ 
tion, Eq. (051): 


■ v4^0 = 


02 / 777 s 2^-3 

|[2(.-l,.i-.b(i) . 

(99a) 


(99b) 

o 

II 

(99c) 


The axisymmetric parts of the non-radiative second- 
order perturbations ((5^^z;“) and ((5^^17) are determined 
by solving the perturbed Euler equation, Eq. (ISTI) . and 
the perturbed mass conservation equation, Eq. (1531) . The 
lowest order in angular velocity contributions to Euler’s 
equation 0 = ((5^^i?a) are given by: 

0 = = -2zuflo(5^>uf) + d^{d^^^U2) 

0= {s'^^E,) = d,{s'^^U2)-iznl(^y'~\ (loob) 

0= = 2tJ7flo((5^^z;r). (100c) 

The integrability conditions for these equations, 
(S^^^Ea) = 0, are given by V[a((5^^E;{,]) = 0. In cylin¬ 
drical coordinates, the lowest-order in angular velocity 
parts of these integrability conditions are 


0 


0 

0 


flnZ /-07 \ 2^-3 


( 

V[,(4)e;^]) = f7o9z(4^0, 

V[^(-54a^]) = floS^ (m{S^^^vT)) ■ 


.Rj 


(101a) 

(101b) 

(101c) 


These conditions, together with the requirement that the 
solution is nonsingular on the rotation axis, determine 

time independent dif¬ 
ferential rotation As before, we denote a par¬ 

ticular choice by 


(4^0 = 0, (102) 

(4>f> = 1)^(1)“"^ (103) 
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The remaining component, is determined from 

the lowest order in angular velocity piece of the perturbed 
mass conservation equation [cf. Eq. (l55|) ]. 

(p(4'^<)) = 0. (104) 

This equation, together with Eq. (I102L shows that the 
only non-singular solution for is 


47rp^ f /ra7\2^-2 

^ 7^ I 4^ 2 (fl) 

-f2flg [ w'5^^'^VL{w')dw' + .(111) 


= 0. (105) 


The scalar parts of the second order non-radiative 
r-mode, (4V) are determined by com¬ 

pleting the solution to the perturbed Euler equation 
( 54 ^ 0 ) = 0, and then solving the perturbed gravita¬ 
tional potential equation. The potential (S^^'^U) is de¬ 
termined by integrating the perturbed Euler Eqs. (llOOal) 
and (|100bl) . Using Eqs. and (11031) we obtain the 
following expression for the axisymmetric part of the so¬ 
lution, to lowest order in angular velocity, 

{2 )rr\ 

= ^\r) 

-f2flg / w'6^^^n{w')dw' + 4^<^ 2 , (106) 

Jo 

where 5^^'^C 2 is a constant. 

The pressure as well as the density perturbations, 5^'^'^p 
and are related to 5^‘^'>U as follows, 


Finally, we use Eq. (IMl) to obtain an explicit formula 
for the second-order differential rotation, 5^^'^VL{w), in 
terms of the second-order radiation-reaction force and the 
second-order velocity perturbation 5)^ v°‘. Of the three 
terms on the right side of that equation, we will see that 
the second and third are higher order in O than the first, 
and we will evaluate the first term to leading order in fl. 

We first use Eq. (155)) to find an explicit form for the 
second-order radiation-reaction force {S^^^/gr}- From 
Eqs. (1551) and (1551) for 4^^^ 

(44g4 = (112 ) 

The second term 4^^'^ is given by Eq. ()103D . 

In the final term, 4^1 ^‘^j by its definition (l50l) . is propor¬ 
tional to a product of components of 4^1^ S^r^v. By 
our initial normalization, = 0 (fl), and we found in 

Sect. HVXI that 4Vis /gr) = 0{l3n'^). 

From Eqs. (1651) . ()112p . and (I103p . we have 


-I--4|y(i)pj(i)p 

P ‘^P 

= 5 ( 2 )$-f p 


P 

' 2 p 2 


7(7 - 2 ) ^ 

p dp 


(J(bp)2^ (107) 


where 7 = d log p/d log p is the adiabatic index. For the 
r-modes, the first-order perturbations and d^^ip are 
O(fl^). So at lowest order in angular velocity, the relation 
between and d^^ip simplifies to 

d(2)[/2 = d(2)$2 + ^d(2)p2. (108) 

p 2 

The gravitational potential is determined by 

solving the perturbed gravitational potential equation, 

y2d(2)^ = qTTd^^V- (109) 


For the r-modes, to lowest order in the angular velocity, 
this equation my be re-written as 

V2d(2)$2 + ^d(2)$2 = ^d(2)[/2. (110) 

IPo IPo 

Using the expression derived in Eq. (11061) for the axisym- 

/o') 

metric part of 5)^112, we find the general equation for 


{5fF<^) 




+(C-i)(i)^ 


(113) 


Equation (11131) implies The second 

term in Eq. (1551) has integrand proportional to (d^d)- 
Because d^d = the integrand is and the 

term itself is O (17^), two orders higher than ( d^^ E . Fi¬ 
nally, the last term in (l64l) is proportional to n(d^^lU“). 
Eq. (ISTI) implies (d^^lU“) = 0(17^), whence the last term 
is again 0(17^). 

With the dominant term in Eq. (1641) determined by 
(d^^E'^), we have 


d 4 17(^7) 


ajzpisfF^) 


(114) 


This integrand can be re-written in a more explicit form 
using Eqs. (11131) and (11031) : 


4^^(^) = (^) 


tx7\ 2^-4 


(£+ 1)2 


(l)‘ 


- 1 ) 


T(to) 


(115) 


2 
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where T{w) is the equation-of-state dependent, mass- 
weighted average of {zjR)'^, 

dz pz^ 

TM = , ■ ( 116 ) 

The limits of integration, ±23(137), in this expression are 
the w dependent values of 2 at the surface of the equilib¬ 
rium star. To lowest order in VL these limits are the same 
as those in a spherical nonrotating star: 

zs(w) = \/ R? — w^. (117) 

The part of the second-order differential rotation that 
is not explicitly caused by the radiation-reaction force, 
(5^puf), is given in Eq. (|103p : 

(!)”-•. (n«, 

Together Eqs. (I115D and (11181) determine (to lowest order 
in n) the time-dependent differential-rotation induced by 
gravitational-radiation reaction: 


2. Equation (11071) determines p 2 k +2 up to the 

( 2 ,) 

ambiguity associated with C 2 k+ 2 - The Pois¬ 
son equation, Eq. 671), with the conditions 
that S)^^^ 2 k +2 vanish at infinity and have no 

monopole part (no change in mass), determines 
(2) (2) 

both S}^ ^ 2 k +2 and the constant S}^'C 2 k+ 2 - 

3. Equation (I107|) (or, alternatively, the Poisson equa- 

f2'l 

tion) gives p 2 k+ 2 , and the equation of state de¬ 
termines 5 ^^'^P2k+2- 

4. Finally, one uses the known first-order perturba¬ 
tion to solve two independent components 

of the curl of the Euler equation, 6 ^^'^ Ea = 0 for 

^N^4k+i and (S^N^k+i) ^las an f{w) 

ambiguity that is resolved by Eq. (IMl) . The fi¬ 
nal component is found from the second- 

order mass-conservation equation. 

C. Secular drift of a fluid element 




„2/3t 


(119) 


The key result of this section is the derivation of an 
explicit expression (11141) for in terms of the 

first-order r-mode. An expression of this kind exists be¬ 
cause the rest of the second-order perturbation, the per¬ 
turbed density, pressure, and potential, are higher-order 
in n. Like the velocity field of the first-order r-mode, 
the second-order differential rotation of the unstable r- 
mode can be found without simultaneously solving for 
the perturbed density and pressure. 

This separation of orders also leads to an iterative 
method for solving the second-order Newtonian pertur¬ 
bation equations at successive orders in fl that mirrors 
the method we have just used to determine the axisym- 
metric parts of at 0 {fl) and 4^’ 4^'!’ 

at 0(11^). At each order, the ambiguity in the Newto¬ 
nian differential rotation is resolved by using Eq. (1641) . 
We assume that the first-order Newtonian perturbation 
equations have been solved to the desired order in fl. We 
suppose one has found the perturbed Newtonian velocity 
4 ^'*^“ to and the scalar quantities in 44 to 

0 (12^^'), and we list the steps to obtain the next-order 

correction: to find 4^^2fe-i-i the scalar quantities to 
f 2l 

1. Because "2fc-i known, and the integrability 
( 2 ) 

conditions = 0 have an additional power 

of n in each term, they are satisfied to at 0(12^^). 
One can then integrate the 177 or 2 component of 
the perturbed Newtonian Euler equation dUl) to 

find 442fe-i-2 up to a constant 442fc+2- 


The differential rotation we have found for the unstable 
r-mode extends the work of Sa and collaborators [i3-[I3 
to obtain the differential rotation of the unstable second- 
order r-mode. The studies of magnetic field wind-up by 
Rezzolla, et al. M, which predated this work, explic¬ 
itly omitted the form of the second order perturbation to 
the velocity field that we have computed here. These au¬ 
thors obtained a secular drift in the position of a fluid 
element by integrating the £ = 2 form of the equations 
for the position <j){t) and 9{t) of a particle whose per¬ 
turbed velocity field is found solely from the first-order 
perturbation 4^'^'* of Eq. dZT]), from the equations 

^ (120a) 

^ (120b) 

The equations are nonlinear in 0(2), ^(2), and the solu¬ 
tion is written to 0{a^). The axisymmetric part of the 
solution is again the part that is not oscillatory in time; 
in our notation, it has the form 


( 0 ( 2 )) = 0 , m)) = a" 




fit. 

( 121 ) 

A secular drift obtained in this way has been used in 
subsequent papers by Cuofano, et al. [nllTl, and by 
Cao, et al. [l^ . 

When one includes the second-order differential rota¬ 
tion of the unstable i = 2 r-mode from Eqs. ()119|) . 
additional terms are added to the secular drift (j){t) of a 
fluid element’s position. The resulting expression is given 
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for t 1 //3 by 




( 122 ) 

Using the expression for in Eq. (I119L with 

Eqs. (I115|) and (11181) . we obtain the following explicit 
form for the second-order drift of an unstable i = 2 r- 
mode: 


im 



1 /t37N 

4 Vfly 


T{w) t. 


(123) 


This expression for the drift {4>{t)) is independent of z, 
and therefore describes a drift that is constant on tjj = 
constant cylinders. The analogous expression for the drift 
found previously by Sa [l^ has this same feature, and 
Chugunov [ 2 ^ observes that the drift in these modes can 
therefore be completely eliminated in the pure Newtonian 
case by appropriately choosing the arbitrary second-order 
angular velocity perturbation. 

For long times (that is, for pt arbitrary but /3 <C U), 
the time dependence t in Eq. (11231) is replaced by (e^^* — 
l)/2/3. This expression is not of order 1//3, but satisfies 
the bound 


e2/3t _ 1 

2/3 


+ 1 

< t -- 


(124) 


for t > 0. 


V. POLYTROPIC STELLAR MODELS 


In this section we evaluate Eq. (11191) , to determine the 
changes in the rotation laws of uniformly rotating poly¬ 
tropes that are induced by the gravitational-radiation 
driven instability in the r-modes. Polytropic stellar mod¬ 
els (polytropes) are stars composed of matter whose 
equation of state has the form 

p = i^pi+i/", (125) 


where K and n, the polytropic index, are constants. We 
start with the simplest case, n = 0, the uniform-density 
models. The only dependence of the differential ro¬ 
tation on the equation of state is in T(-n7), the 

mass-weighted average of (z/R)^ at fixed w defined in 
Eq. (jll6[) . This average can be evaluated analytically in 
the uniform-density case: 


T(to) 


— Tjp' 
3i?2 


4(^) 

3i?2 


(126) 


Combining this result with Eqs. (I115L (11181) and (I119L 
gives 




n 


/tZ 7 \ 2^-4 


{e + m-7) 
24 

(£ 2 - 1 ) / ^ 



(127) 



FIG. 1: Differential rotation from the £ = 2 r-mode 

instability evaluated on a cross section through the rotation 
axis of a slowly-rotating uniform-density star. The solution 
scales with time as 


In particular, for the 1 = 2 r-mode, the radiation-reaction 
induced differential rotation has the form 


(5(2)0 = O 


3 

2 



5 

8 




(128) 


which is positive in a neighborhood of the poles and 
negative near the equatorial plane. Figure [T] illustrates 
the gravitational-radiation driven differential rotation 
i5(2)0/fl from the i = 2 r-mode instability of a slowly- 
rotating uniform-density star. This figure shows contours 
of constant ^( 2 ) 0 / 0 , on a cross section of the star that 
passes through the rotation axis. For example, this figure 
ilustrates that (5^2)U/U k, —9/8 near the surface of the 
star at the equator. This indicates that the angular veloc¬ 
ity of the star is reduced by an amount ~ — (9/8)fla2e2(^‘ 
in this region, where is the amplitude of the r-mode, 
and U is the angular velocity of the unperturbed star. 
Similarly this figure illustrates that i5(2)U/U k, I near 
the poles. The angular velocity of the star is enhanced 
by the r-mode instability in these regions. 

The equilibrium structures of n = 1 polytropes can 
also be expressed in terms of simple analytical functions, 
but the integrals that determine T(tx7) in Eq. (11161) can 
not. We therefore evaluate these quantities for all the 
n 0 polytropes numerically. 

The structures of the non-rotating Newtonian poly¬ 
tropes are determined by the Lane-Emden equations, 
which are generally written in the form, 

I 

where 9 is related to the density by p = Pc9'^, with 0=1 
at the center of the star and 0 = 0 at its surface. The 


6 
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variable ^ is the scaled radial coordinate, r = a^, with 

( 130 ) 

47rG 

We solve Eq. (I129D numerically to determine the Lane- 
Emden functions O(^), use them to evaluate the density 
profiles of these stars, p(r) = PcQ^, and finally perform 
the integrals numerically in Eq. (jll6l) that determine the 
mass weighted average T(n7) of (z/R)"^ for spherical poly¬ 
tropes. Eigure [5] illustrates the results for a range of 
polytropic indices. Because they are more centrally con¬ 
densed, stars with softer equations of state, i.e. poly¬ 
tropes with larger values of n, have smaller T(ci7). This 
is most pronounced near the rotation axis of the star 
where to = 0 and values of z'^ in the dense core dominate 
the average. Figure [3] illustrates from Eq. (I115j) . 

the differential rotation induced by the gravitational- 
radiation driven instability in the i = 2 r-modes for 
polytropes having a range of polytropic indices n. This 
graph shows that the equatorial surface value (to = R) 
of is the same for all the poly tropes. This is 

not a surprise, because T(to) = 0 there for all equations 
of state. Stars composed of fluid having stiffer equa¬ 
tions of state, i.e. smaller values of n, have larger values 
of near the rotation axis where to = 0. Fig¬ 

ure |4] illustrates the differential rotation induced by the 
gravitational-radiation induced instability in the r-modes 
of n = 1 polytropes having a range of different spherical 
harmonic mode index I values. The figure portrays a 

/ n \ 

differential rotation (5 jy “/fl induced by gravitational ra¬ 
diation that, like the magnitude of the linear mode, is 
more narrowly confined to the equatorial region near the 
surface of the star as the r-mode harmonic index i is 
increased. 

VI. DISCUSSION 

The radiation-reaction force uniquely determines the 
exponentially growing differential rotation of the unsta¬ 
ble, nonlinear r-mode. We have found expressions for 
the rotation law and for the corresponding secular drift 
of a fluid element and have obtained their explicit forms 
for slowly rotating polytropes. The formalism presented 
here describes an r-mode, driven by gravitational radi¬ 
ation reaction, at second order in its amplitude a, and 
restricted to a perfect-fluid Newtonian model. We now 
comment briefly on the meaning of the work within a 
broader physical context. 

First, a realistic evolution involves coupling to other 
modes, because realistic initial data has small, nonzero 
initial amplitudes for all modes and, at higher orders in 
a, other modes are excited by the r-mode itself. As a 
result of the couplings, the r-mode amplitude will sat¬ 
urate, and studies of its nonlinear evolution (see [1, @ 
and references therein) suggest a saturation amplitude of 
order 10“"^ or smaller. By the time the mode reaches sat¬ 
uration, the amplitude of daughter modes may be large 



FIG. 2: Dimensionless ratio of the integrals T(to) defined 
in Eq. (fTTfil) that determines the gravitational-radiation in¬ 
duced differential rotation in polytropic stellar models having 
a range of polytropic indices n. 



FIG. 3: Differential rotation induced by the gravitational- 
radiation instability in the 1 = 2 r-modes for a range of poly¬ 
tropic indices n. 


enough that their own second-order axisymmetric parts 
contribute significantly to the differential rotation law. 

Second, when there is a background magnetic field, the 
growing axisymmetric magnetic field generated by the r- 
mode’s secular drift can change the profile of the growing 
differential rotation [2^. The second-order Euler equation 
(1451) is altered by the second-order Lorentz force per unit 
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FIG. 4: Differential rotation induced by the gravitational- 
radiation instability in various r-modes of n = 1 polytropes 
for a range of spherical harmonic mode index I values. 


mass, given in an ideal magnetohydrodynamics approxi¬ 
mation by ((5(2) /magnetic) = Q;^ (V X i?) X B]). 
This will be of order the radiation-reaction force after an 
amplitude-independent time^ 


t ' 




lO^’s 


P 10®G 


lO^^g/cm^ Bq 


(131) 


where Ia is the Alfven time associated with the back¬ 
ground held, Ia = y/Airp/BQ. After this time and until 
the mode reaches its nonlinear saturation amplitude, we 
expect that the radiation-reaction force will continue to 
drive growing differential rotation. The functional form 
of this differential rotation, however, will be determined 

by both (5(2) and ((5(2) /magnetic)- 


^ For a magnetic field that grows linearly in time, we have 

/magnetic) 

AirpH 

The second-order radiation reaction force is given by 
~ a^pflR, implying that the Lorentz force 
q ,2 { 5 ( 2 ) has comparable magnitude after a time given 

in Eq. JTsTJ. Here we follow Chugunov [ 2 ^ . Chugunov uses 
this argument to conclude that the magnetic field will not be 
significantly enhanced after it reaches B ~ 10®(a/10“^)^ G, 
but his analysis is restricted to the case where the gravitational 
radiation-reaction force on the r-mode is negligible. We have 
checked the conclusion of continued growth for Shapiro’s model 
of a uniform-density cylinder with an initial magnetic field 1271 . 
by adding a forcing term of the form of the second-order axisym- 
metric radiation-reaction force |2^ . We expect the amplification 
factor of the magnetic field to be limited by the value of the mode 
amplitude, at nonlinear saturation, not by the value of the 

field, unless the initial magnetic field is of order 10^^ G or larger. 


After nonlinear saturation, we expect the growth of dif¬ 
ferential rotation and of the magnetic field to stop within 
a time on the order of the Alfven time. This is because 
( 1 ) the radiation-reaction force is now time independent, 
and ( 2 ), with a background magnetic field, there should 
no longer be a zero-frequency subspace of modes associ¬ 
ated with adding differential rotation. Reason (2) means 
that the differential rotation and the magnetic field at 
the time of mode saturation become initial data for a 
set of modes whose frequencies are of order the Alven 
frequency. The second-order axisymmetric part of the 
r-mode after saturation becomes effectively a system of 
stable oscillators driven by a constant force. Such sys¬ 
tems have no growing modes, and therefore no secularly 
growing magnetic field. 

The explicit form of the secular drift we obtain is new, 
but its mag nitude is consistent with that used in earlier 
work ( 8 I-I 12 I [isjl that examines the damping of the un¬ 
stable r-mode by this energy transfer mechanism. This 
damping mechanism becomes important whenever the 
rate of energy transfer to the magnetic field (by winding 
up magnetic field lines or, for a superconducting region in 
a neutron-star interior, by stretching magnetic-flux tubes 
or other mechanisms), is comparable to the growth rate 
of the unstable r-mode. Assuming the energy transferred 
to the magnetic field is not returned to the r-mode and 
that a large fraction of the core is a type II supercon¬ 
ductor, Rezzolla et al. Q estimate that the instability 
will be magnetically damped for a magnetic field of order 
10^2 Q_ As noted above, we expect this magnetic damp¬ 
ing mechanism to play a role only if the magnetic field 
reaches this 10^2 Q range prior to nonlinear saturation 
of the r-mode. We think it likely that a limit on mag¬ 
netic field growth imposed by saturation means that this 
field strength can be reached only if the initial field is 
not far below 10^2 Q. In addition, for an initial field of 
order B > 10^2 Q or larger, if all axisymmetric pertur¬ 
bations that wind up the magnetic field have frequency 
higher than or of order the Alfven frequency, we conjec¬ 
ture (based on the toy model mentioned in Foonote 4) 
that the enhancement of the magnetic field will be too 
small to damp the r-mode. 

Finally, if the magnetic field is large enough to signifi¬ 
cantly modify the structure of the first order r-modes, all 
of the calculations here would need to be modified. Previ¬ 
ous studies, however (l^. [29l - l33| . find that field strength 
B > 10^"^ — 10^® G is needed to significantly alter the 
linear r-mode of a star with spin greater than 300 Hz. 
When the viscous damping time is comparable to the 
gravitational-wave growth time, one would also need to 
include viscosity in the 2 nd-order equations that deter¬ 
mine the differential rotation. 
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The symbols in Table U are listed by order of appear¬ 
ance in the paper, starting with Sec. We omit a few 
symbols that are used only where they are defined. 


Symbol 

Q 

p 

V 

p 

$ 

h 

U 

E°- 

fan 

j£m g£m 

Ni 

Q 


^( 1 ) x( 2 ) 

On ? On 




UJn 

tpN 

P'^b 


P 

(SQ) 


g(2) 

o \TC 


Q 




{S^n^F) 

Mo, R 

flo 

n 

Qn 

fan 


TABLE I: Notation 
Meaning 

the set of variables {p, v, p, <!>} 
mass density 
fluid velocity 
fluid pressure 

Newtonian gravitational potential 
fluid specific enthalpy 
effective potential 

E°' = 0 is the Newtonian Euler equation 

radiation-reaction force 

mass and current multipoles 

a constant defined in Eq. m 

fluid angular velocity 

rotational symmetry vector xy — yx 

dimensionless amplitude of r-mode 

first-order perturbation of Q: d^gl^^Q 

second-order perturbation of Q: 1 daQ\^_Q 

first- and second-order Newtonian perturbations 

(no radiation reaction) 

■uj, z dependence of perturbation: Eqs. (|3^-p5ll 
a correction in first-order perturbation due to 
radiation reaction 

subscript ± denotes even (+) or odd ( —) t/i-parity 
under the diffeomorphism (j> ^ 211 — (j) 
frequency of Newtonian r-mode 
ijjN = WNt + mcj) 

projection operator orthogonal to (f>: Eq. 
imaginary part of frequency of unstable r-mode 
axisymmetric part of 5Q 
2nd-order Newtonian perturbation with a 
particular choice of 5),^ Q{-w) 
arbitrary function of zu in second-order 
Newtonian differential rotation 
second-order differential rotation, 
radiative part of second-order perturbation 
W defined in Eqs. and m 

effective driving force for (S^^v): Eq. (1571) 
mass and radius of spherical stellar model 

VW/iP 

dimensionless angular velocity, fl/flo 
part of g that is nth order in Q: Eq. dzoi 
part of fcR orthogonal to 
function for which 5^^'^ fcR — 


Appendix B: Gravitational Wave Energy and 
Angular Momentum Fluxes 

The expression for the radiation reaction force foR 
given in Eq. o was derived by constructing a force that 
reproduces the standard expressions for the time aver¬ 
aged gravitational wave energy and angular momentum 
fluxes: 



The expression given here for the angular momentum 
flux, Eq. (IB2I) . is somewhat more compact than the 
standard post-Newtonian expression (cf. Thorne [l^ 
Eq. 4.23). We express this flux in terms of the mag¬ 
netic type mass and current multipole moments, 
and S'J™, which we define as 

/j"* = NeVITT J (B3) 

pr^{v-rxV)Ys*^^d^x. (B4) 

These magnetic type mass and current mutipole mo¬ 
ments can be expressed in terms of the standard 7^"* 
and 

f* £m 


Q £rn 




- - ^/(£-m)(£ + m+l)I^ ^+^(x + iy) 

— — iy) 


(B5) 




i_ / (£ + m)(£-m+l) _ ... 

2 \/ £{£+l) ^ ^ 


zJ Z, 


^/W+T) 


(B6) 


where x, y and z are unit vectors. Both of these ex¬ 
pressions are based on the following identity for vector 
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spherical harmonics: 


y “*' Im 
B 


2 \/ l{l + l) 


i / (£ + m)(£-m+l) ^^^ 
2 \/ i{l+l) 




^/W+^) 


-.Y^ 




(B7) 


Using this transformation, Eq. (IB2I) reproduces the 
standard post-Newtonian expression (cf. Thorne [l^ 
Eq. 4.23). The calculation needed to verify that the 
expression for the radiation reaction force fcR given in 
Eq. d?]) satisfies the time averaged gravitational wave en¬ 
ergy and angular momentum flux expressions given in 
Eqs. (IBll) and (IB2I) is straightforward, but lengthy. 


Appendix C: Integrating 5p 

For rotating equilibrium stellar models having poly¬ 
tropic equations of state with polytropic index n, the 
density p oc (distance to the surface)" near the star’s sur¬ 
face. We assume here that the surface of the perturbed 
star is smooth as a function of a and x. Although the 
surface itself is smooth, the behavior of p near p = 0 im¬ 
plies that VaP diverges for n < 1 and S/aS^bP diverges for 
n < 2. It follows that 5^^^p and p diverge because they 
involve first and second derivatives, respectively, of the 
unpertubed density. We show, however, that continuity 
and differentiability of the star’s surface as a function of 
X and a imply finiteness of the integrals pdz and 

1^00 pdz, when 5^'^^p is regarded as a distribution. 

We first verify the claimed behavior of p for the unper¬ 
turbed polytrope and then use the form of the Lagrangian 
perturbation of the enthalpy to deduce the behavior of 
p and p near the surface. Denote by zf (a, t, w, cj)) 
the values of z at the top and bottom parts of the sur¬ 
face of the perturbed star. We again introduce the poly¬ 
tropic function 0 , related to the specific enthalpy by 
9 = po/[{n + l)po\ h. Then p = po9^Q{zg - z)Q{z - Zg), 
where 6 ( 2:5 — z) = 1 for Z 5 > z and 0 (z 5 — z) = 0 
for Z 5 < z. For the unperturbed rotating polytrope, 9 is 
finite with finite derivatives at the surface of the star.® 
The lack of smoothness in p at the surface thus arises 
from the fact that n is not an integer. We now show for 
the perturbed polytrope that p is again proportional to 
(distance to the surface)" to second order in a. 


^ For the unperturbed star, Eq. (Il3]l implies 9 = po/{{n+l)po\(€ — 
■!> + where £ is the constant injection energy per unit 

mass. Caffarelli and Friedman prove that p is Holder continuous, 
p g C®’“(R®) [ 3 ^, and the Poisson equation then implies G 
C'^’“(R®). Thus 9 has one-sided first and second derivatives at 
the surface. 


The vanishing of 9 at the surface of the perturbed star 
is equivalent to the vanishing of the Lagrangian pertur¬ 
bation of 9 at the unperturbed surface: 

A0 = 0, (Cl) 

where 

A0 := 0 ( 0 , t, if-I-^) — 0(0, t, x), (C2) 

with ^(a, t, x) the exact Lagrangian displacement - a vec¬ 
tor from the position x of each fluid element in the un¬ 
perturbed star to its position if-|-^ in the perturbed fluid. 
Our assumption that the surface changes smoothly as a 
function of a and x is then the requirement that ^ and its 
derivatives are smooth at the surface of the star. Writing 

{= +0(a^) (C3) 

and taking derivatives of (IC2I) with respect to a, we have 

5 ( 1)0 = A(i)0-^(i)“Va0, (C4a) 

j(2)0 ^ A(2)0-^(2)“Va0-C(^)“Va5(i)0 

_ l^(i)a^(i)fcv,Vfc0. (C4b) 

Then 5(^)0 and 5(^)0 are hnite at the unperturbed surface, 
and, to second order in a, we can write for 0 the Taylor 
expansion 

0(a, z, to) = dz9\^+ (z — Z 5 ) -I- 0(z — Zg}'^, (C5) 

for z < Zg. The corresponding expansion for p = po9”' is 
thus 

p(a,z,TO) = po{-dz9\^+)'^{z^-z)^ + 0{zj-z)'^+^. 

(C6) 

We can now show directly that the integrals 
J^^S^^'^pdz and J^^S^^^pdz are finite for polytropic 
equations of state with any polytropic index n > 0 for 
which the equilibrium star has a finite surface. More 
precisely, they are finite everywhere except the equator, 
where the range of integration vanishes. 

For a given value of to, we choose with 0 < Z'^ < 
Zg and 0 > Z~ > Zg for all a < e, for some hnite e > 0. 
We write the integral as a sum of three parts, 

/ OO nOO pZ 

Spdz= / 6pdz+ 5pdz+ Spdz. (C7) 
-00 J Z~ J Z'^ J — 00 

In the hrst integral on the right side, p and 5(^)p are 
hnite, so we need only consider pdz, S^^^pdz, 

and the corresponding integrals near the bottom part 
of the surface. Because the hniteness argument is 
identical for the integrals near Zg and Z 5 , we consider 
the integrals near Z 5 . 
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We have 

daP = da^[po{-d;,e\^+)'^{zj - z)'^ 

+ Oiz+-zr+^]Q{z+-z)} 

= Po(-9^6»|^+)”aa[(4 - 2:)”0(4 - Z)] 

+o{zt - zr 

= -pod^zU-d.9^rd4iz+ - zreiz+ - z)] 
+0{zt - zr, (C8) 

implying 

s^^^p = -Poi^^^^r-drrUsdzK^s - zTQizs - z)] 

+0{zs - zr, (C9) 

where we have used the relation daZg\a=o = 

From Eq. (Oil . we have 

dip = -po5a4(-^^^L+)"'9«9z[(4 - - ^)] 

+0(4 - z)"-i 

= prdo^ztfi-dArTdliizt - z)"0(4 - ^)] 

+0(4-z)"-\ (CIO) 

implying 


5 ( 54 Q)- The perturbed Euler equation, Eq. (l95|) . with 
radiation-reaction force then has the form 

= [(wAr+£fl)5a6 + 2iflVa<('b]5^^^U^ 

= zv, (4^4 - 4^4). (Di) 

Using the slow-rotation form (ISTl) of a ;at and Eq. (13) of 
Ref. [ 3 ^ , we write the solution to this equation for 
as 

4^^;“ = (44 - 44 ) , (D2) 

where the inverse of Q~^ is the tensor Q°'^ = 
with 

= ~WW) + l)"V“zV'>z - *(£ + 1)V“4] . 

43 ) 

With 6 ^v°‘ replaced by the expression on the right side 
of Eq. (ID2I) . the mass conservation equation, Eq. (IMl) 
becomes an elliptic equation for 44 — namely 


Va 


pg“4b (44-44)' 

+^^^4^(44 - 44) = *u/344. (d4) 


s^^^p = Ipoii^^^^rH-drrhsdiiizs - ^)" 0(4 - ^)] 

+ Oizs-zr-^. (Cll) 

Finiteness of f S^^dz is immediate from the integra- 
bility of (zs — zr~^ for n > 0. For had to retain 

the factor 0 (zs — z), and we kept it for as well to 
display pairs of analogous equations. From Eqs. (IC9I) 
and (|C11I) . the leading term in each of and d^^^ is 
a ^-derivative, and we immediately obtain the integrals 

poo 

/ d^^'>pdz = p{Z)r^'>^U,+0{zs-Z+r+\{CU) 

Jz+ 


r d^^rdz = !^^^^ir^)%r"+oizs-z+r. 

Jz+ z Zs — z^ 

(C13) 


The integrals jr^d^^dz and jr^d^'^rdz are therefore 
finite as claimed. 


Appendix D: Ordering in Q, of 44 

To make the heuristic argument of Sect. IIV Al more 
precise, we use the two-potential formalism of Ipser and 
Lindblom [s^ to write an explicit form for in terms 

of 44 and 44 . Because that formalism uses the 
complex version of a perturbation, we write 44 = 


The potentials 44 and 44 are determined by this 
equation, together with the Poisson equation, 

= 47rp4(4)[/ _ 4 )$)^ (D5) 

and the two boundary conditions, 


lim = 0 

r^oo ^ 

(D6) 

and 


= (44|s - 44 + I44bh) 

^ = 0; (D7) 

here S is the surface of the unperturbed star and the 
Lagrangian displacement 4^ is defined by 

4 ^ • 

(D8) 


Using the value of a; at from Eq. m, and Eq. (lD2]) , we 
can write the second boundary condition as 


Wb(4) u - 44) + (4^ u - 4 ) $) = 0 . 

(D9) 

To find the orders in U of 44, 44 and 44, we be¬ 
gin with the relations 44 = C)(U^) and, from Eq. (|M|) . 
4^4 = 0{p.f5). From the Poisson equation (ID5I) . 44 
and 44 are the same order in U. From Eq. (ID4I) . 
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we then have = 0{n'^d‘'l^T) + 0{Vt^P) = 

0{n^l3). Then 

<5^$ = C>(4^V) = OiS’l^T) = 0{np). (DIO) 
Finally, Eq. (ID2I) implies 

= C>(D-id 3/3) = C>(D2/3). (Dll) 
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